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^Department of Physics and Astronomy, University of Glasgow, 
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We present new results on T and J 1^ spectroscopy using Lattice NRQCD. Charmonium results at quenched /? 
= 5.7 show agreement with experiment within understood systematic errors. We compare bottomonium results at 
quenched (is 6.0 and 5.7 and show that the IP— 15 splitting scales with Ay, provided a correction is applied for the 
0{a?) errors inherent in the gluon field configurations. There is a clear scale difference between charmonium and 
bottomonium results in the quenched approximation. We compare quenched results at /9=6.0 for bottomonium 
with those obtained on HEMCGC configurations at /9=5.6 using 2 light flavours of dynamical quarks. We show 
that extrapolations to k/ =3 are consistent with experiment for the ratio of the 25 — 15 splitting to that of the 
IP — 15. We extract a value for the T-r]b splitting extrapolated to ra/ =3. 



1. Introduction 

The NRQCD collaboration has been active for 
a few years in the field of heavy quark physics. 
Our aim is to test QCD on the lattice by compar- 
ing charmonium and bottomonium spectroscopy 
to experiment as accurately as we can. This has 
required : 

• The development of a non-relativistic effec- 
tive action appropriate to the physics [1]. 

• A perturbative improvement scheme to re- 
move discretisation errors order by order 
[2,3]. 

• High statistics calculations with multi- 
exponential fits to multiple smeared corre- 
lators [4-6]. 

Heavy-heavy spectroscopy lends itself very well 
to these techniques. An important feature is that 
sensible estimates can be made of the remaining 
systematic errors. 

Previous results have focussed on T spec- 
troscopy on quenched lattices at /3=6.0 from the 
Staggered collaboration [6]. The new results 
reported below are on quenched lattices at [3 
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= 5.7 (from UKQCD) and on lattices with 2 
light fiavours at /3=5.6 (from HEMCGC). View- 
ing the NRQCD action as an expansion in pow- 
ers of I we have used, as before, leading 
and next-to-leading spin-independent terms along 
with leading spin-dependent terms, v'^/c^ k, 0.1 
for hh and 0.3 for cc. Discretisation errors in the 
leading order terms are also removed. The coeffi- 
cients of all terms have been set to their tree-level 
values once the link fields have been 'tadpole- 
improved' by dividing by the fourth root of the 
plaquette [2,3]. 

2. Charmonium results 

The charmonium spectrum as determined on 
quenched configurations at /3=5.7 is shown in Fig- 
ure 1 [7]. The bare lattice quark mass is fixed (at 
0.8 in lattice units) from the requirement that the 
mass in the non-relativistic dispersion relation for 
the rjc should be correct. The expected system- 
atic error is 30 MeV from v^' j c^' terms, and this 
error is clearly evident in the fine structure of the 
spectrum. Quenching corrections could also ap- 
pear at this level in the hyperfine splitting. 

The inverse lattice spacing, determined from 
the spin-averaged \P — \S split- 
ting is 1.23(4) GeV. By spin-averaged splitting we 
mean Pi) - 0.25 x \im{^Si) ^ m{^S^)\. The 
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2S — IS splitting gives a consistent value but is 
not well determined because excited states decay 
very quickly at low /3. 

3. Upsilon results 

The bottomonium spectra that we obtain on 
quenched configurations at /3=6.0 and «/=2 con- 
figs at /3=5.6 (HEMCGC) are shown in ref. [8]. 
The expected systematic error here is only 5 MeV. 
There are clear signs that the 25—15' and 1_P— 15 
splittings are inconsistent in the quenched ap- 
proximation. This is a refiection of the scale- 
dependence of the inverse lattice spacing because 
the quenched coupling constant runs incorrectly. 
This problem is much improved on the «/=2 con- 
figurations. 

We can also compare the IP — 15 splitting ob- 
tained at different values of /3. Figure 2 shows 
this splitting divided by Ay obtained from the 
plaquette [9] vs Ay in lattice units. For asymp- 
totic scaling the results at /3 =6.0 and 5.7 should 
lie on a horizontal line. They do this, but only 
after a correction has been applied for the effect 
of 0(0^) effects from the plaquette gluon action. 
This correction can be calculated perturbatively 
[8] or non-perturbatively by comparing solutions 
in corrected and uncorrected heavy quark poten- 
tials [10]. This figure also makes clear that differ- 
ent values are obtained from different quanti- 
ties, even when discretisation errors are corrected. 
a-i(T) > a-i(*) > a-Hjrip). 

Returning to a comparison of quenched and un- 
quenched results, we show in Figure 2 the vari- 
ation with of the ratio of the 25 — 15 split- 
ting to that of the IP — 15. All splittings are 
corrected perturbatively for the effect of gluonic 
0(0^) errors. The corrections amount to only 2% 
at /3=6.0. It is clear that the ratio is ri/ depen- 
dent. A linear extrapolation takes us to the ex- 
perimental value of 1.28 at the value of nf ap- 
propriate to the momenta inside an T of 3. This 
is very gratifying and shows clearly that sensible 
values for physical quantities (e.g. Ug [8]) can 
be obtained by extrapolation in nf. It would be 
useful to check the results for nj=A. 

In a similar way the hyperfine splitting m(T) — 
m(rii,) can be extrapolated to nf =3. We obtain a 
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Figure 1. NRQCD simulation results for the char- 
monium spectrum plotted relative to the spin- 
average of the J/^ and the ryc(15) using an in- 
verse lattice spacing of 1.23 GeV from the 1_P— 15 
splitting. Experimental values are indicated by 
dashed lines. Error bars are shown where visible, 
and only indicate statistical uncertainties. 



value of 44 MeV. The expected systematic error 
of 5 MeV dominates the error in this result. It 
may be true that «/=4 is more appropriate to 
this short distance quantity and then one should 
correct it upwards perturbatively. 
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Figure 2. Asymptotic scaling of the IP— IS split- 
ting for hh for quenched results in the V scheme. 
Open circles represent raw results at /3 = 5.7 and 
6.0 and filled circles represent results corrected for 
gluonic 0(0^) errors. Filled circles are offset for 
clarity and the horizontal line is drawn to guide 
the eye. Two other points appear on the graph - 
the square from the cc spectrum at /3 = 5.7 and 
the diamond from the a ^ limit of nip from 
ref. 11. These two points have been rescaled by 
the ratio of the physical values of 1_P — 15 for cc 
or nip to the IP — IS splitting of hh (taken to be 
452 MeV). 
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